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Abstract 

Bifurcation with symmetry is considered in the case of an isotropy subgroup 
with a two-dimensional fixed point subspace and non-zero quadratic terms. In gen- 
eral, there are one or three branches of solutions, and five qualitatively different 
phase portraits, provided that two non-degeneracy conditions are satisfied. Condi- 
tions are also derived to determine which of the five possible phase portraits occurs, 
given the coefficients of the quadratic terms. The results are applied to the problem 
of bifurcation with spherical symmetry, where there are six irreducible representa- 
tions for which the subspace of solutions with cubic symmetry is two-dimensional. 
In each case, the number of solutions and their stability is found. 

1 Introduction 

The topic of bifurcation with symmetry has wide-ranging applications, including buckling 
of rods, convection patterns in fluids and the structure of tumours and embryos. In the 
physical context, a highly symmetric state becomes unstable as a control parameter is 
varied, leading to a state with reduced symmetry. If there is a high degree of symmetry, 
then the stability problem is degenerate, in the sense that several eigenvalues pass through 
zero simultaneously as the parameter is varied; the number of eigenvalues is equal to the 
dimension of one of the irreducible representations of the original symmetry group. This 
means that the dynamics in the neighbourhood of the bifurcation is governed by a number 
of coupled, nonlinear equations which in general cannot be solved. 

This complicated multi-dimensional problem can be simplified by restricting attention 
to the flow-invariant fixed-point subspace of a subgroup of the original symmetry group. 
If this subspace is one-dimensional, then the problem reduces to a single differential 
equation, so that the usual methods of bifurcation theory can be applied, leading to a 
bifurcation of transcritical or pitchfork type. This is essentially the equivariant branching 
lemma (Vanderbauwhede 1980, Cicogna 1981, Golubitsky, Stewart and Schaeffer 1988): 
if the fixed-point subspace is one-dimensional, then a unique branch of solutions exists in 
the vicinity of the bifurcation. For example, in the problem of the buckling of a rod with 
square cross-section, solutions exist in which the rod buckles either parallel to two of the 
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sides or diagonally. The equivariant branching lemma is known to generalize to the case 
where the dimension of the fixed-point subspace is odd, but it is in general untrue when 
the dimension is even. 

In this paper, the case of a two-dimensional fixed-point subspace is considered, following 
on from earlier work. Mcleod and Sattinger (1973) studied this problem, and showed 
that there are one or three solution branches for generic values of the coefficients of the 
quadratic terms, corresponding to five different phase portraits; however they regarded 
the double zero eigenvalue as a degeneracy rather than a consequence of symmetry. 
Lauterbach (1992), Leblanc et al. (1994) and Lari-Lavassani et al. (1994) formulated the 
problem in the context of equivariant bifurcation theory and obtained specific bifurcation 
theorems. In section 2 below, the basic theory of the subject is recalled, together with a 
summary of some of the known bifurcation theorems. Section 3 gives specific conditions 
under which bifurcation can be guaranteed, in terms of the coefficients of the quadratic 
terms, and also gives the relationship between these coefficients and the five possible 
phase portraits. In section 4 these results are applied to the problem of bifurcation from 
spherical symmetry, 0(3), in the fixed point subspace of the symmetry group of the cube, 
which is two-dimensional for six of the irreducible representations of 0(3). 



2 Branching in two dimensions 

There is no simple generalization of the equivariant branching lemma (EBL) to the case 
of a two-dimensional fixed point subspace. In this section, some of the key definitions of 
the subject of bifurcation with symmetry are recalled and some of the known results are 
summarised. 

Near a multiple bifurcation point, the eigenfunctions form the basis of a vector space V. 
The symmetry group G acts on V by multiplying points in V by matrices that form a 
representation of G. In most cases this representation is irreducible, meaning that there 
is no proper G-invariant subspace of V . The symmetry of a point v in V is described by 
its isotropy subgroup, defined by {g G G : gv = v}. 

For any subgroup H C G, the fixed-point subspace of H is 

Fix(if) = {v G V : hv = v for all h G H}. (1) 

Fixed-point subspaces are important because they are invariant subspaces: if an initial 
condition lies in Fix(iJ) then the system remains in Fix(iJ) for all time. 

The normalizer N(H) is defined by 

N(H) = {geG: g-'Hg = H}. (2) 

This is the largest subgroup of G in which H is normal. Equivalently, this is the subgroup 
of G which maps Fix(iJ) to itself: 

N(H) = {geG:gye Fix(H) for all y G Fix(if)}. (3) 

Clearly H C N(H), and if N(H) = H then there is no symmetry within Fix(H), since 
all elements of H act as the identity in Fix(H). In general, the symmetry group acting 
on Fix(iJ) is the quotient group N(H)/H. 
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Henceforth it will be assumed that the bifurcation problem with a symmetry group G 
satisfies the following conditions. 



1. G is a compact Lie group acting absolutely irreducibly on a vector space V. This 
implies that Fix(G) = or Fix(G) = V. The latter case only holds for the trivial 
one-dimensional representation of G, so we may assume that Fix(G) = 0. 

2. The function / : R x V — > V is smooth and G-equivariant, i.e. f(X,gv) = gf(X,v) 
for all g 6 G, v 6 V. This implies that /(A, 0) = gf(X,0), and hence /(A, 0) G 
Fix(G), so /(A, 0) = 0. Since the action of G is absolutely irreducible, the lineariza- 
tion of / is a scalar multiple of the identity, h(X)I. 

3. h(0) = 0, and h'(0) ^ 0, i.e. there is a stationary bifurcation at A = and the 
eigenvalue passes transversely through zero as A passes through zero. 

4. H is an isotropy subgroup of G with Dim(Fix(if )) = 2. 

5. The quadratic terms in the Taylor expansion of / do not all vanish in Fix(H). 

Note that the first three conditions are the standard ones for the application of the 
EBL (Ihrig and Golubitsky 1984, Golubitsky, Stewart and Schaeffer 1988, Chossat and 
Lauterbach 2000). 

The equations /(A, v) = within Fix(if) can then be scaled, for small A, by A = e/i//i'(0), 
v = ex, to give, as e — > 0, 



involving six constants a, b, c, p, q, r; from the last of the above conditions, not all of 
these are zero. 

These two equations were discussed by Mcleod and Sattinger (1973), who showed that 
for generic values of the six constants there exist either one or three solutions. However, 
a general bifurcation theorem cannot be deduced from this, because in the context of 
equivariant bifurcation theory the six constants are not independent. There are two 
types of constraints on these coefficients (see, for example, Stewart and Dias 2000). 
The first type of constraint arises from normaliser symmetries which act on Fix(if) as 
shown in (J3J); for example, there may be a normaliser symmetry acting as a reflection 
y — > —y, in which case p = b = r = 0. Secondly, there may be other constraints arising 
from equivariance with respect to the original symmetry group G; these are sometimes 
referred to as 'hidden' symmetries. 

An example to illustrate this point is the symmetry group of the tetrahedron including 
reflections, which is isomorphic to the symmetric group S4 and also the group O of 
rotations of the cube. In the natural three-dimensional irreducible representation, the 
group is generated by permutations of the three Cartesian coordinate axes and a sign 
change of any two, and the bifurcation equations consistent with these symmetries to 
quadratic order are 











fix + ax 2 + bxy + cy 2 , 
\iy + px 2 + qxy + ry 2 , 



(4) 
(5) 



= [ix + yz = fxy + xz = jjz + xy. 
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Note that there is only one equivariant quadratic term; this is a feature common to many 
bifurcation problems, including the cases of Sn symmetry and 0(3) symmetry, although 
there are also examples where this is not the case, such as icosahedral symmetry (Hoyle 
2003). The fixed point subspace of the symmetry group Z 2 generated by y <-> z is 
two-dimensional and the bifurcation equations reduce to 

= /ix + y 2 = ixy + xy. (6) 

The normaliser is Z 2 x Z 2 , including the element (y, z) — > (—y, —z), and so the normaliser 
symmetries force p = 5 = r = 0in(|nj), but we also have a = and c = q, arising from 
the original symmetry group. Furthermore, there are only two non-zero solutions to (JHJ), 
x = —[/,, y = ±fi, not one or three as expected from (JJJ |3J). (In fact, in this example, a 
third solution in which y — z — does exist when cubic terms are included, as required 
by the EBL since this subspace is one-dimensional.) 

Some rigorous results on bifurcation in two-dimensional subspaces are as follows. A result 
of Lauterbach (1992) (Theorem 1.2), obtained using degree theory, is that at least one 
non-zero solution exists if, in addition to the five conditions above, 

g- 1 (0) = 0, i.e. Q = 0^ x = y = 0, (7) 

where Q represents the quadratic terms in (JU|3J). 

Lauterbach's condition (JZj) is sufficient but not necessary. Leblanc et al. (1994) showed 
(Theorem 2.1) that in the case of gradient dynamics, when Q can be written as the 
gradient of a scalar quantity, a branch of solutions to the truncated problem (jjj [3J) always 
exists. Furthermore (Theorem 3.2), they showed that a smooth branch of solutions of 
the non-truncated problem exists if either <5 _1 (0) 7^ or a fourth-order polynomial in 
the six coefficients a ... r is non-zero. 

Hence, by combining the results of Lauterbach (1992) and Leblanc et al. (1994), bifur- 
cation to non-zero solutions always occurs occurs in the case of gradient dynamics, since 
both the cases Q _1 (0) = and Q _1 (0) 7^ are covered. 

3 Phase portraits 

In this section, solutions to the system (jH EJ) are considered in more detail. First, the 
condition under which there are one or three roots is derived. 

Multiplying (jlj) by y, J5J by x and subtracting yields the cubic equation 

= cz 3 + (b - r)z 2 + (a - q)z -p (8) 

for z = y/x, which generically has one or three solutions. The condition for the existence 
of three roots depends on the quantity 

P = 27cV - 18cp(b -r)(q-a)- Ap(b - rf - Ac(q - a) 3 - (q - af{b - r) 2 . (9) 

There are three roots if P < 0, one root if P > 0, and a double root if P = 0. An 
alternative derivation of the existence of one or three solutions, employing the implicit 
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function theorem, is given by Mcleod and Sattinger (1973). Given a solution for z = y/x, 
substitution into (@J) gives either the trivial solution x = or x = —fi/(a + bz + cz 2 ). 
This fails to give a finite solution if the denominator is zero, and it is straightforward to 
show that this occurs if the quantity 

R = (ar — pc) 2 + (cq — br)(aq — bp) (10) 

is zero. R is the resultant of the two quadratic polynomials appearing in (jU EJ. It is a 
well known result of algebraic geometry that these two homogeneous polynomials have 
a non-zero solution if and only if the resultant is zero (see, for example, Cox, Little and 
O'Shea, 1998). Thus the non- degeneracy condition R ^ is exactly the condition (JJJ), 
Q _1 (0) = 0. The degenerate case R — corresponds to a solution at infinity. 

The non-degeneracy conditions giving one or three solutions are therefore P / and 
R 7^ 0. These conditions have not been given explicitly in earlier work, although Leblanc 
et al. (1994) give P in the special case of gradient dynamics. 

It is also of interest to consider the two-dimensional phase portraits within Fix(H), by 
considering the time dependent form of (JUG)}: 

x = fix + ax 2 + bxy + cy 2 , (11) 
y = fiy + px 2 + qxy + ry 2 . (12) 

In this system, the straight line passing through the origin and any fixed point is an 
invariant line, since if z — y/x satisfies (JBJ) then z = (yx — xy)/x 2 which is zero, since 
this is proportional to the same cubic appearing in (|S|). This observation forces all fixed 
points to be saddles or nodes, and prohibits the existence of any periodic orbits. The 
invariance of these lines is a consequence of the truncation, and is broken at higher order. 

Stability of the fixed points is determined by the eigenvalues of the Jacobian 

j_/fi + 2ax + by bx + 2cy \ , 
V 2px + qy fi + 2ry + qx J ' 

One eigenvalue is — fi, corresponding to the eigenvector lying along the invariant line. 

In general, there are five possible different types of phase portrait (McLeod and Sattinger 
1973). Which of these occurs is determined by the values of P and R and a third quantity, 

I = aq + br — bp — cq. (14) 

The following result gives the conditions on P, R and / that determine the phase portrait. 

Theorem 

In the system (fmH2|l . 

(a) If P > and R > there is one fixed point, which is a saddle point. 

(b) If P > and R < there is one fixed point, which is a node. 

(c) If P < and R > there are three fixed points, one node and two saddles. 

(d) If P < 0, R < and I > there are three fixed points, two nodes and one saddle. 

(e) If P < 0, R < and I < there are three fixed points, all of which are saddles. 

Proof 

Note first that P, R and I are invariant under rotation of coordinate axes. This can be 
verified by determining that under an infinitesimal rotation through an angle 9, 

da db dc 

Te =b + P > de = q + 2c - 2a > de =r - b ' 
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de 11 de ' r d# 

To check that I is invariant, we expand 
— = q(b + p) + (r — p) (g + 2c — 2a) — q(r — 6) — b(q — a) + (a — c) (2r — 2p — b) — b(c + q) 

de 

and verify that this quantity is zero, and similarly for P and R. 

Exploiting the rotation invariance and the fact that there is always at least one non-zero 
fixed point for P ^ 0, R ^ 0, we can always rotate coordinates so that a fixed point lies 
on the x axis with x < 0. After this rotation, it follows from (fT2*|) that p = and from 
(fTTj) that the fixed point is at x = —fi/a. It can be assumed that \i > 0, because the case 
[A < is equivalent under a sign change of x, y and t, and so a > 0. The eigenvalues of 
the fixed point at x = —fi/a, y = are, from pfljl . — /x and — g/a), so the fixed point 
is a saddle point if a > q and a node if a < q. Consider now the five cases above in turn. 

(a) If P > and R > 0, then after rotating so that p = and a > 0, P > gives 

4c(a - g) > (6 - rf (15) 
so that either c > and a > g, or c < and a < q. The condition i? > becomes 

ar 2 + cq 2 — brq > 0, (16) 

which can be written as 

(a — q)r 2 + cq 2 — rg(6 — r) > 0. 

Now if c < and a < q, this quadratic function of r must have real roots for the inequality 
to be satisfied, so (6 — r) 2 > 4 (a — g)c, which contradicts (fT5J). Thus the only possibility 
is c > and a > g, so the fixed point is a saddle point. The phase portrait in this case 
is shown in figure 1(a). 

(b) If P > and R < then (fTK|) holds and R < can be written as 

(a — g)r 2 + cq 2 — rq(b — r) < 0. 

Using the same argument as in the previous case, if a > q and c > the inequality can 
only be satisfied if (b — r) 2 > 4(a — q)c, contradicting ()15j) . Therefore c < and a < q, 
and the fixed point is a node, as shown in figure 1(b). 

(c) If P < and R > there are three fixed points. Using the result of Mcleod and 
Sattinger (1973) that at least one of these is a saddle point, a rotation of coordinates can 
be carried out so that p = and a saddle point lies at x — —fi/a, y = 0; hence a > q. 
The quadratic equations for the coordinates (xi, yi) and (x 2 , y 2 ) at the other fixed points, 
obtained from (JH) and (JHJ), are 

(ar 2 + eg 2 — brq)x 2 + (r 2 — 6r + 2cq)p,x + c/j 2 = 0, (17) 
(ar 2 + eg 2 — brq)y 2 + (2ar — rg — bq)[x,y + (a — g);U 2 = 0. (18) 

Since i? > 0, the leading term in each quadratic is positive, and so from (fTHj) . ?/i?/2 > 0. 
Thus all three fixed points lie in the same half plane, y > or y < 0. In the case that all 



lie in the upper half plane, consider an anticlockwise rotation that brings the next fixed 
point onto the negative axis. After this rotation, the two points not on the x axis have 
Z/1Z/2 < 0, so a < q and so the point on the axis must be a node. After a further rotation 
so that the third fixed point lies on the x axis, both other fixed points have y < 0, so 
a > q and the point on the axis is a saddle. Hence for P < and R > 0, there are three 
fixed point in the same half plane, and the middle of the three is a node while the other 
two are saddle points (figure 1(c)). 

(d) If P < 0, R < and I > 0, then 4c(a - q) < (b — r) 2 and a 2 r 2 + (cq - br)aq < 0, so 
cq — br and aq must have opposite signs. Since / > 0, aq + br — cq > and so both aq 
and br — cq must be positive. Since P < 0, there are three fixed points. After a rotation 
of coordinates, one lies on the x axis and the other two satisfy (fT7|ITR|) . At the other two 
fixed points, one eigenvalue is — \i and the other, from (fT3j) . is 3/i + (2a + q)x + (2r + b)y. 
The sum of the 'non-radial' eigenvalues at these fixed points can be found using the sums 
of the roots obtained from (|17llT%j) . The result is qafi(r 2 + b 2 — 2br + Acq — 4ac) / R, which 
is negative. Thus, at least one of these fixed points must be a stable node. Given that 
one fixed point is a node, this point can be rotated onto the x axis and an argument 
exactly analogous to that of case 3 can be followed, leading to the conclusions that the 
three fixed points lie in the same half plane, the middle one is a saddle point and the 
other two are nodes (figure 1(d)). 

(e) If P < 0, R < and / < then, as in case (d), cq — br and aq must have opposite 
signs. In this case aq + br — cq < 0, so aq < and br — cq < 0. Hence a 2 — aq > and 
so the fixed point on the x axis is a saddle point. Thus in this case all three fixed points 
are saddles, as shown in figure 1(e). 

Remarks 

1. The five possible phase portraits are illustrated in figure 1, for fi > 0. The case 
fi < is obtained by reversing all the arrows and rotating through 7r. 

2. Stable fixed points can only occur in cases (b), (c) and (d). Of course, these 
stability assignments only apply within Fix(H), and there are directions transverse 
to Fix(if ) which may be stable or unstable. In fact, it is generally the case that all 
solution branches obtained are unstable in the full system, since it can be shown 
that one eigenvalue is — fi but the sum of the eigenvalues is rifi, where n is the 
dimension of the representation (see, for example, Chossat et al. 1990). 

3. As P passes through zero there is a saddle-node bifurcation, leading to a transition 
from case (a) to case (c) or from case (b) to case (d). As R passes through zero, 
a fixed point passes through infinity, giving transitions between cases (a) and (b), 
(c) and (d) or (c) and (e). There can be no direct transition between cases (d) and 
(e) because / = is incompatible with R < 0. 

4. In the case of gradient-like dynamics, there are additional constraints b = 2p and 
q = 2c. Under these conditions it can easily be shown that P > => R > and 
R < =^> / < 0, so that there are only three possible phase portraits, (a), (c) and 
(e). 
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(a) (b) (c) 




(d) (e) 

Figure 1: The five possible types of phase portrait for (fTT|IT2|) . when fi > (for \i < 0, all 
arrows are reversed). The open circle indicates the trivial solution and filled circles are 
non-zero solutions. The five case are (a) one solution, saddle, (b) one solution, node, (c) 
three solutions, saddle, node, saddle, (d) three solutions, node, saddle, node, (e) three 
solutions, all saddles. 
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4 Application to bifurcation with spherical symme- 
try 



An interesting example of bifurcation with symmetry is the problem of bifurcation from 
an initial state with spherical symmetry. This problem has many physical applications 
and has been studied at length by Sattinger (1978, 1979), Ihrig and Golubitsky (1984), 
Chossat et al. (1990), and Matthews (2003). The irreducible representations have dimen- 
sion 21 + 1 and the eigenf unctions are the spherical harmonics of degree I. Of particular 
interest are patterns with the symmetry groups of rotations of the Platonic solids: the 
tetrahedron, T, with 12 elements, the cube / octahedron, O, with 24 elements, containing 
T, and the icosahedron / dodecahedron, /, with 60 elements, also containing T. When 
/ is odd, the quotient group N(H)/H always includes the group Z% , generated by point 
inversion through the origin, corresponding to a rotation through n if Dim(Fix(if)) = 2, 
so the quadratic terms are all forced to be zero by the normaliser symmetry. For even 
/, the isotropy subgroups are all of the form H © Z 2 , where H is a subgroup of SO(3) 
and © indicates a direct product. The groups O © Z% and / © Z% are maximal isotropy 
subgroups, and satisfy the condition N(H) = H . The dimensions of their fixed-point 
subspaces are [1/4] + [1/3] -1/2 + I for O © Z 2 C and [1/5] + [1/3] -1/2 + 1 for I © Z 2 C , 
where [x] is the largest integer less than or equal to x. 

The dynamics is gradient-like, so the following two existence results follow directly 
from the observation at the end of section 2, by determining the values of I for which 
Dim(Fix(#)) = 2. 

1. There exists at least one, but not more than three branches of solutions with 
isotropy O © Z 2 C for I = 12, 16, 18, 20, 22, 26. 

2. There exists at least one, but not more than three branches of solutions with 
isotropy I © Z 2 C for / = 30, 36, 40, 42, 46, 48, 50, 52, 54, 56, 58, 62, 64, 68, 
74. 

Since the dynamics is gradient- like, the phase portraits within the two-dimensional sub- 
space Fix(H) may have any of the three forms shown in figure 1(a), (c) or (e) provided 
that the non- degeneracy conditions P ^ 0, R ^ are satisfied. For the case I = 12, 
H = O © Z% ', the computation of the coefficients was carried out by Leblanc et al. 
(1994); the result is, after removing some common factors, 

mi3V&m _ 7350i v / 8398 _ l^A^llm ssgiigv/ssis 

390 ' °~ 4940 ' P ~ 3 ' 7 ~~ 4370 ' 

with b = 2p and q = 2c. Calculating the values of P, R and I shows that all three are 
negative. Hence there are three distinct branches of solutions with isotropy O © Z% for 
/ = 12, and within this subspace, all three are saddle points. 

A similar calculation was carried out for the other values of / for which Dim(Fix(0 © 
Z2)) = 2, using the computer algebra package Maple to do the cumbersome manipula- 
tions. First, the 21 + 1 equations for the amplitudes of the spherical harmonics Y^iO, <p), 
m — —I . . .1 are constructed using the fact that the quadratic coefficients are the Clebsch- 
Gordan coefficients (Sattinger, 1978). The subspace Fix(D 4 © Z$) of dimension 1 + [1/4] 
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is found simply by restricting to those modes in which m is a multiple of four. The sub- 
space Fix(0 © Z2 ) is then obtained from this by imposing invariance under a rotation of 
it/2. The results of these calculations are as follows. For I = 12, 16, 22 and 26, P, R and 
/ are all negative, so there are three stationary solutions and each solution is a saddle 
point in Fix(0 © Z%), as in figure 1(e). For I — 18, P < and R > 0, so there are three 
solutions, one of which is a node and the other two are saddles (figure 1(c)). For I = 20, 
P > and R > 0, so there is only one solution, which is a saddle point. 
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